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Abstract. We study dilepton production in proton-proton, Cu + Cu as well as in Au + Au collisions at 
the center-of-mass energy ^/s= 200 GeV per participating nucleon pair within an extended statistical 
hadronization model. In extension to earlier studies we incorporate transport calculations for an estimate 
of uncorrelated e^e~ -pairs from semileptonic D meson decays. While the invariant mass spectrum of 
dielectrons is well understood in the p +p collisions, severe discrepancies among different model scenarios 
based on hadronic degrees of freedom and recent data from the PHENIX Collaboration are found in 
heavy- ion collisions in the low mass region from 0.15 to 0.6 GeV as well as in the intermediate mass 
regime from 1.1 to 3 GeV when employing the standard dilepton sources. We investigate, furthermore, 
the background from correlated dileptons that are not emitted as a pair from a parent hadron but 
emerge from semileptonic decays of two correlated daughter hadrons. Our calculations suggest a sizeable 
contribution of such sources in central heavy-ion collisions in the low mass region. However, even the upper 
limits of our calculations are found to be far below the dilepton mass spectra of the PHENIX Collaboration. 
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1 Introduction 

The goal of this vifork is to study the production of corre- 
lated electron-positron pairs in proton-proton and heavy- 
ion collisions at -ysj^f^=200 GeV. We will evaluate the 
invariant mass spectrum diV/dmg+e- within an extended 
hadronization model up to invariant masses of 4 GeV thus 
covering the low mass as well as the intermediate mass and 
charmonium regime. We confront our calculations with 
corresponding measurements of the PHENIX collabora- 
tion (within the experimental acceptance) which has taken 
data [1,2] up to 4 GeV in TOg+e- at mid-rapidity. In the 
case of proton-proton collisions the PHENIX collabora- 
tion has found out that the measured spectrum can be 
described very well up to masses of 4 GeV, if one takes 
into account all relevant hadronic sources of dileptons in 
the analysis. This was done by a simultaneous measure- 
ment of all hadron rapidity densities and transverse mo- 
mentum spectra around mid-rapidity; by using these ex- 
perimental rapidity densities one can then estimate the 
dilepton yields at different invariant masses due to the 
known hadronic decays to e^e~ pairs. Independently, mi- 
croscopic transport calculations within the Hadron-String- 
Dynamics (HSD) framework [3,4] have come to the same 



conclusion when incorporating the measured cross section 
for cc pairs from the PHENIX collaboration [5] . 

In contrast to the proton-proton collisions, the mea- 
sured invariant mass spectrum of dileptons in Au -\- Au 
collisions have so far not been properly understood the- 
oretically. Instead, theoretical estimates are found to de- 
viate up to a factor of 4 or 5 from the PHENIX data for 
central Au-|-Au collisions in the low mass regime [5,?] (cf. 
also Fig. 42 in Ref. [2]) and by up to a factor 2 to 3 for in- 
termediate masses. Such discrepancies among (hadronic) 
models and experimental data have not been observed at 
lower Super-Proton-Synchrotron (SPS) beam energies and 
in different collision systems [5,7] where a major broaden- 
ing of the vector- meson resonances is reported. This might 
suggest that non-hadronic dilepton channels could be re- 
sponsible for the discrepancies observed so far. This issue 
needs further independent investigations. 

In this work, we will partly repeat the analysis by the 
PHENIX collaboration in the intermediate mass range es- 
pecially with respect to the contribution from charmed 
meson decays but instead of using the measured yields 
as input for the dilepton sources, we wish to implement 
different models in order to calculate the rapidity den- 
sities of different hadrons and subsequently estimate the 
emission of dileptons from their decays. These models are 
controlled by the PHENIX data for p -)- p collisions. 
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We recall that the invariant mass spectrum of dilep- 
tons can be divided in three fairly distinct regions in each 
of which different physics processes are dominant. Below 
the (/) meson mass, the region referred hereafter as the 
low mass region (LMR: mf,+ e- € [0.0 ; 1.1] GeV), the 
dilepton production is dominated by the decays of non- 
charmed mesons, i.e. mesons with essentially light quark 
content (M,d, s). In the intermediate mass region (IMR: 
nif^+f^- g[1.1 ; 3.2] GeV), i.e. in between the (j) meson 
and J/'f/' mass, the invariant mass spectrum of electron- 
positron pairs is dominated by the semileptonic decay 
products of open charm mesons. Strictly speaking this 
is just a background for the "true" dilepton sources, but 
since this component is always present in the analysis, 
one needs to carefully evaluate the contribution from open 
charm as well. Furthermore, above about 3 GeV of invari- 
ant mass the direct decays of charmonia become dom- 
inant and provide a constraint on the number of pro- 
duced cc-pairs - forming bound states - once the char- 
monium suppression is controlled independently. We con- 
centrate in this article on studying the LMR and IMR 
thus exceeding previous approaches that focused on the 
LMR and extrapolated to the IMR. The high mass region 
(HMRiTTie+e- >3.8 GcV) of the dilepton invariant mass 
spectrum is dominated by the Drell-Yan process and B 
meson decays which we will not address here. 

In describing the yields or ratios of particle yields of 
hadrons consisting of u, d and s quarks phenomenologi- 
cal models, in particular "thermal models" , have proven 
to be very useful due to their simplicity and low number 
of adjustable parameters. In this work, we will evaluate 
the yields of light mesons within the statistical hadroni- 
zation model which has been applied to high-energy ele- 
mentary [8-10] and especially heavy-ion [11-21] collision 
experiments in order to calculate the yields of different 
hadron species with fairly a lot of success. We mention 
that the CERES Collaboration has also analyzed their 
dilepton data in the LMR on the basis of the statistical 
hadronization model with some success at Super-Proton- 
Synchrotron (SPS) energies for Ph + Au collisions [22]. 

Unlike for the bulk meson production in the LMR, the 
statistical hadronization model can not be applied to esti- 
mate the yields of charmed hadrons in the IMR. Instead, 
we need to rely on experimental information here in order 
to estimate the differential yields of charmed hadrons. To 
this aim we will formulate an 'extended statistical model' 
including early charm-pair production, collective flow of 
the hot and dense matter as well as correlated and un- 
correlated semileptonic decays of D mesons. The amount 
of D meson rescattering will be followed up by the HSD 
transport approach [23] in order to estimate the amount of 
surviving correlated semileptonic decays in the heavy-ion 
collisions. 

In extension to earlier studies we will, furthermore, 
explore the contribution of correlated semileptonic decays 
stemming from kaon-antikaon pairs that emerge from the 
decay of heavy parent hadrons. The most important of 
these sources are proportional to 7^, where 75 denotes 
the strangeness suppression factor which is low in proton- 



proton reactions but close to unity in nucleus-nucleus col- 
lisions. Our upper limits for these contribution will be 
compared to the PHENIX data within the proper accep- 
tance windows as well as the conventional dilepton sources 
mentioned above. 



2 An extended statistical hadronization model 

The statistical hadronization model (SHM) has been ap- 
plied successfully in calculating the number of emitted 
hadrons in high-energy collision systems [8-21] . This model 
is well documented in the references given above and ac- 
cordingly we will introduce the main concepts only. We 
evaluate the hadron yields in the grand-canonical ensem- 
ble because the calculations simplify substantially if one 
does not require exact conservation of Abelian charges 
and/or energy-momentum. Our choice is motivated by 
two reasons: First of all experimental observations show 
that the grand-canonical ensemble is sufficient enough, i.e. 
that at RHIC energies the data are well described under 
the approximations we have chosen. The other reason is 
that in order to evaluate the hadron yields in the canon- 
ical ensembles, one needs to know the volume as well as 
the exact (integer) charges on an event by event basis. 
However, the PHENIX collaboration has measured only a 
small fraction of the emitted hadrons while a large part 
of the system is never observed. We would need to make 
severe assumptions for the part of the system not mea- 
sured, if we were to implement the canonical formalism 
for the calculation. We also note that the canonical effects 
arc most pronounced for heavy and exotic hadrons, while 
the bulk of the dilepton emission arises from the low-mass 
mesons, which are produced abundantly and do not suffer 
from canonical suppression effects. Thus, we deem that 
performing the analysis in the grand-canonical ensemble 
should be good enough for our purposes and indeed this 
seems to be confirmed by our work (see below). 

In the SHM, the primary hadron multiplicity of hadron 
type i is calculated (in the on-shcll Boltzmann approxima- 
tion) according to 



N,,^V 



2 J, -hi 



l^n. ^l^-^./T ^-^p2+m^/T ^3^^ (^) 



In Eq. (1) Jj denotes the spin, p the momentum and rrii 
the mass of the particle while qi is a vector consisting of 
the baryon, electric and strangeness charges of the hadron 
species i. The state of the "thermal" fireball is specified by 
its temperature T, volume V and chemical potentials (col- 
lected in the vector fi) for baryon, electric and strangeness 
charges. 

Several independent experimental measurements have 
verified that the mid-rapidity region is actually almost 
net charge free in proton-proton and Au + Au collisions at 
-y/ij,ji^=:200 GeV [24-26]. The baryon chemical potential 
is expected to be of the order of 30 MeV [19] in central 
Au + Au collisions and we have used this value throughout 
at all centralities while the fxs and ij,q are set to zero. 



Manninen et al.: Dilcpton production in p + p, Cu + Cu and Au + Au collisions at ^/s = 200 GeV 



Once the chemical potentials are fixed, there are three 
free parameters characterizing our system: the tempera- 
ture, the strangeness under-saturation parameter 75 and 
the overall normalization volume V . The auxiliary param- 
eter 7^ is necessary to include in our analysis in order to 
take into account the empirical fact that the strange par- 
ticle yields are strongly suppressed with respect to SHM 
estimates in elementary particle collisions. We have chosen 
75=0.6 in our analysis in accordance with statistical had- 
ronization model fits [25,27] to proton-proton collisions at 
this beam energy. In Au -\- Au collisions, the 7s parameter 
increases monotonically from the p-\-p value to unity as a 
function of centrality [19] and the effect will be discussed 
in detail in the forthcoming sections. The common nor- 
malization factor V for all hadron species is determined 
in different collision systems from the measurements (see 
below) while for the temperature we have chosen the value 
T = 170 MeV in all systems based on the SHM fits to 
proton-proton and Au + Au data at this beam energy. 

We have included the same collection of hadron species 
in our analysis as has been included in the works quoted 
above from which we have taken the thermal parameters. 
The mean primary hadron and resonance yields of each of 
the hadron species included in the analysis are calculated 
according to Eq. (1). For resonances with width larger 
than 2 MeV, Eq. (1) is convoluted with the relativistic 
Brcit-Wigner distribution and integration over the mass 
and momentum is enforced. Once the mean primary yields 
are known, we assume that, from event-to-event, the mul- 
tiplicity distribution of each species is governed by the 
Poisson distribution characterized by the mean multiplic- 
ities evaluated with Eq. (1). We then sample the Poisson 
distributions for each of the hadron species in each event to 
obtain the primary multiplicities and choose momenta for 
every hadron according to the Boltzmann distribution. All 
unstable resonances are then allowed to decay according to 
the most recent branching fractions taken from the parti- 
cle data group tables [28] . This way we know the momenta 
of every final state particle and thus it is straightforward 
to take into account any geometrical and kinematical ex- 
perimental cuts. 

We note here that the statistical hadronization model 
is useful only for evaluating the relative yields of different 
hadron species since the rapidity and transverse momen- 
tum spectra of all hadrons emitted in the high energy col- 
lision experiments do not resemble thermal distributions. 
This is not a problem for us as long as our results do not 
depend explicitly on the details of the spectra. Indeed, 
this is the case for the "true" dilepton emission from a 
single decaying hadron, i.e. the invariant mass of the lep- 
ton pair does not depend on the momentum of the parent 
hadron and neither does the (Lorentz invariant) number of 
hadrons. Thus, we may evaluate the number of produced 
hadrons and dileptons within the statistical hadronization 
model even though the spectra are not correct. 

Unfortunately, the discussion above holds true only if 
the measurement is performed in 47r, i.e. if all hadrons 
are measured or at least if the experiment extrapolates 
the measured kinematical region to the unmeasured re- 



gions as well. The PHENIX dilepton measurement is car- 
ried out in a narrow rapidity window around y = Q and 
only leptons with px > 200 MeV are taken into account. 
We have also taken into account the azimuthal geometry 
and effects of the magnetic field on the charged leptons 
of the PHENIX detector [2] , namely, that electrons and 
positrons (of charge q) are accepted in case both of the 
conditions 

(t^min < (t) + q < 0max ; fec = 0.206 rad GeV/c 

Pt 

0min < </» + q < 0max ] kmcH = 0.206 rad GeV/c 

Pt 

are simultaneously fulfilled. The PHENIX detector con- 
sists of two arms with the angular coverage 0niin 
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11 
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16 



Ig^: 'f'max — IQ^ ' 



The limited acceptance in rapidity is not a severe prob- 
lem here due to the approximate boost invariance of the 
systems around mid-rapidity at top RHIC energies. In 
order to make the rapidity distributions of non-charmed 
hadrons wider, we have randomly boosted^ (event by event) 
our " fireball" along the beam axis so that the rapidity dis- 
tributions of pions become compatible with the BRAHMS 
measurements [29]. 

On the other hand, the limited acceptance in px raises 
some problems, because the statistical hadronization model 
tends to over-populate the low pT part of the spectrum 
compared with the experimental distributions and thus 
too few leptons hit the PHENIX acceptance window of 
Pt > 200 MeV. We have solved this problem by assuming 
that the created clusters' transverse momentum is nor- 
mally distributed (with mean /ip^=0 but (pt) > 0) and 
fitted the width of the clusters' pT distribution together 
with the system volume V to the PHENIX data [30,31] in 
p + p collisions and in 11 different centrality classes in the 
case of Au + Au collisions. Comparison of the data and 
model calculations are shown in Figure 1 for the p+p (top 
panel) and central Au + Au collisions (bottom panel) . The 
description of the data is similar at all centralities. The 
resulting widths for the clusters' transverse momentum 
distributions are shown in Figure 2 while the scaling vol- 
umes (divided by the volume in proton-proton collisions) 
are shown in Figure 3. 

The dilepton yields are measured in wider centrality 
classes than the pt spectra from which we have deter- 
mined the (Tpj, and V and thus in the following calcula- 
tions for dielectrons we have used interpolated values for 
ap^ and V based on the values shown in Figures 2 and 3. 



3 Decay widths 

In the low invariant mass region the dominant sources of 
(correlated) dileptons are the direct and Dalitz decays of 
light mesons. The dielectron decay channels taken into 
account in this analysis are listed in the Table 1. Each of 



^ assuming a Gaussian distribution with width ay — 4.2 
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Fig. 1. (Color on-line) Transverse momentum spectrum of vr , 
7r~, K^ , K~ , p and p [30,31] in proton-proton collisions (top 
panel) and (0-5%) most central Au + Au collisions (bottom 
panel) at y^j^j^=200 GeV compared with our model calcula- 
tions. The width of the clusters' transverse flow profile along 
with the system volume are fitted to the data shown in the 
figure. Kaons have been divided by a factor of 10 and nucleons 
by a factor of 100 for clarity. 



the direct decays results in a sharp peak in the mass spec- 
trum at the meson nominal mass while the Dalitz decays 
yield a continuum spectrum from zero invariant mass up 
to the mass of the decaying meson. Let us note here that 
there are many other hadrons and their resonances decay- 
ing radiatively into dileptons than the ones listed in Table 
1. These could be (and are) important in different kinds 
of collision systems. For example in heavy-ion collisions 
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Fig. 2. (Color on-line) The width of the clusters' transverse 
momentum distribution in 11 centrality classes in Au + Au col- 
lisions at -yS[^j.j=200 GeV (filled spherical symbols) as a func- 
tion of the number of participants. The square symbol denotes 
the width in p + p collisions. 
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Fig. 3. (Color on-line) The system volume normalized to the 
volume in p + p collisions in 1 1 centrality classes in Au + Au 
collisions at -ysj^j^=200 GeV as a function of the number of 
participants. Deviations from the A'^part scaling (straight line) 
are visible only in the central collisions. 



at low beam energies, the Dalitz decay of A resonances 
dominate the low mass region of the dilepton invariant 
mass spectrum (see e.g. [32-34]). Above mid SPS beam 
energies, however, the number of mesons exceed the num- 
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ber of baryons in heavy-ion collisions and at RHIC beam 
energies, the emission of dileptons from baryons is over- 
whelmed by orders of magnitude by the mesonic sources. 
Thus, we do not consider the dileptons stemming from de- 
cays of baryons in this work because the contribution is 
completely negligible at all invariant masses. 



The decay widths for the direct decays of vector mesons 
depend on the mass of the decaying resonance. In prac- 
tice this matters for p° direct decays only, since all other 
mesons are sufficiently narrow that the decay width can 
be considered constant. For the p" direct decay, the decay 
width reads [36] 
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V 


ip' — >• e+e" 


rp' — > 7e"'^e" 


- 


D mes. 


- 


- 


D* -^o^iye + X 



Table 1. List of decay channels relevant for dielectron pro- 
duction in p-|-p collisions at y^=200 GeV. For D mesons, 7 
semileptonic (electron(positron) + anti-neutrino(neutrino) + 
1 or 2 light hadrons) decay channels are considered. For the 
neutral (D" and D^) mesons and D^, there are 6 semileptonic 
decay channels taken into account 



The decay probability of a meson into a pair of lep- 
tons depends on the invariant mass of the lepton pair. A 
generic expression for the decay probability is known from 
Ref. [35] 
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Here mi , mx and M are the masses of the lepton, the de- 
caying meson and the invariant mass of the dilcpton pair, 
respectively. The form factors F^{M) have been studied 
extensively both experimentally and within different mod- 
els. In this work we have employed the form factors arising 
from the vector-meson dominance model considerations [4, 
35] 



f')^77(7V/) 



f 



M' 



(0.72GcV)2 



pIT 



'(M) = 1 + 



5.5 



rM^ 



Thew- 
dr'^ 



GeV2 
7r"/+Z" channel is calculated from 

r 



>7r" i+r 



''^ ''2a 



dM 



M 



X 



1 



M2 



4to2m2 



4mf 
M2" 

3/2 



1 



2mf 
M2" 



with the form factor 
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in which M is the p -meson mass and niQ denotes the pole 
mass. The mass dependent branching fraction of the p'^ 
meson into a pair of leptons is obtained from Eq. (6) by 
dividing it with the mass dependent total width of the p'^ 
meson: 
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(7) 



Above the (f> meson mass, the dilcpton invariant mass 
spectrum attains contributions mainly from the decays of 
charmed hadrons. Radiative decays of J/ip 's into dilep- 
tons have been studied in detail in [37]. Since it is not pos- 
sible to (completely) disentangle the direct {J/tp -^ l^l^) 
and Dalitz decays {J/ip — > 7 Z+Z") of J/ip in the collision 
experiments, one needs to take into account also the Dalitz 
decay of 3/ip in the analysis. This will modify somewhat 
the shape of the invariant mass spectrum of the dilep- 
tons stemming from decays of charmonia. We have imple- 
mented the analytical formula [37] 
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(8) 



which has been used successfully in describing the spectral 
shape of the radiative decays of 3 ftp measured both in 
DESY as well as in PHENIX p+p collisions [1]. In Eq. (8) 
mx is the mass of the decaying particle, M the invariant 
mass of the dilcpton pair and r = -y/l — Ami/M'^. This 
distribution diverges when E^ — )■ (M — )■ mx) and thus 
a cut in energy must be introduced. A suitable value for 
the mass cut-off has been found [37] to be around ii^min ~ 
10 MeV, which we have also employed. Integrating Eq. 
(8) gives us the widths of the radiative charmonia decays: 

p. I / -4,-^^ c+c- ^ Q 32 r-^/'/'^e+e- g^j^j ^.^-'^7 6+6- ^ 034 

F"^ ~^° ° . We have used these widths in evaluating the 
branching ratios for the charmonia Dalitz decays in our 
calculations. It is worth mentioning that our estimated 
branching fraction BR(J/-0^- 7e+e~) is about twice the 
value listed in the most recent PDG book [28]. However, 
our choice agrees somewhat better with the shape of the 
dielectron spectrum near the 3/ip peak than in the case of 
the PDG value for the 3/^p Dalitz decay branching ratio. 
The branching fraction for the ip' Dalitz decay is not yet 
measured and thus a comparison is not possible. 

For all other hadrons - not explicitly mentioned above 
- we have assumed a relativistic Breit-Wigner spectral 
function and the partial widths are then evaluated in a 
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simplified procedure taking into account only trivial mass 
threshold effects for the different decay channels to correct 
for the available phase space. 



4 Charmonium and continuum background of 
dileptons from heavy quark decays 

The weak decays of D mesons D ^i- I + i/i + h, in which h 
denotes one or two non-charmed hadrons, constitute the 
main source of the " dilcpton continuum" in the interme- 
diate mass region at RHIC energies. At RHIC energies 
in proton-proton collisions, there is most often zero or a 
single charmed quark-anti-quark pair created. When this 
pair of charmed quarks hadronizes, the most likely result 
is that each of the charmed quarks end up in a I? and 
D meson. When the D mesons subsequently decay into 
leptons and hadrons, we may have an extra lepton pair 
stemming from the D meson decays in the final state. It 
is very difficult to subtract the leptons originating from 
the D meson decays in the collision experiments and thus 
the measured dilepton invariant mass spectrum usually in- 
cludes this contribution. The contribution of the D meson 
decays is only significant far away from the true sources of 
dileptons and dominates the spectral shape between the 
'peaks' of the (j) and 3/tp mesons. Thus, one needs to care- 
fully consider the D mesons as dilepton emitting sources 
in high-energy collision experiments. 

The slope of the dilepton continuum in the IMR arises 
as a superposition of the momentum distributions of the 
measured leptons coming from the individual D meson de- 
cays. Since the measurement is carried out at mid-rapidity 
and we are interested in invariant masses larger than 1 
GeV, it is clear that the major contribution to the in- 
variant mass of the dilepton continuum in the IMR arises 
from the transverse momenta of the decaying open charm 
mesons. Accordingly, it is important to model the trans- 
verse momentum spectrum of the D mesons more care- 
fully than in the longitudinal direction, i.e. the rapidity 
distribution. 

We will employ a very simple model to evaluate the 
rapidity distributions of charmed hadrons. Namely, we 
assume that all charmed quark-anti-quark pairs are pro- 
duced via splitting of a hard gluon created in the initial 
hard collisions. We also assume that the hard gluons - from 
which all cc-pairs originate - are created via gluon-gluon 
fusion processes in the initial hard scatterings of the glu- 
ons from the target and projectile. In this case, the final 
charmed hadron rapidity distributions will closely follow 
the rapidity distribution of the hard gluons emitted in the 
collision experiment. 

In order to fix the rapidity distributions of charmed 
hadrons we harness the idea of limiting fragmentation [38] , 
which has been verified experimentally at ultra-relativistic 
beam energies both in hadron- hadron [39,40] as well as in 
heavy-ion collisions [41]. Let us consider a collision of two 
gluons with momentum fractions of xi and X2 of the collid- 
ing projectile and target. One can show (see e.g. [42]) that 
the resulting parton rapidity distribution at large momen- 
tum fraction (i.e. xi ^ X2 or vice versa) is proportional 



to the parton distribution function itself 
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and is approximately independent on the Q^ scale due 
to Bjorken scaling. In Eq. (9) pr and y are the trans- 
verse momentum and rapidity of the produced gluon, ttin 
is the mass of the beam particle (in this case the pro- 
ton). Thus, in order to estimate the rapidity distribu- 
tion of the charmed hadrons, we have to adopt a suitable 
parametrization for the parton distribution function in Eq. 
(9). We have chosen the following NNLO pQCD best fit 
parametrization from Rcf. [43] for our gluon distribution 

xg{x) ^ x°-{l- x)''{l + ji^/x + j2x) with 
a = -0.118 b = 9.6 71 = -3.83 and 72 == 8.4. (10) 

This parametrization is given at Q^=9 GeV^ « -^j/i/. ^^"^ 
we approximate the gluon distributions at higher Q^ with 
the same parametrization. 

We can now evaluate the rapidity distribution of charmed 
hadrons with a modified Brodsky-Gunion-Kuhn (BGK) 
model [44] introduced in Rcf. [45], in which the parton 
number density of produced partons along the beam axis 
is proportional to a triangle defined by the momentum 
fractions xi and X2 as follows: The center of mass of the 
colliding partons move with rapidity Tjcm = atanh( ^^~i^^ ). 
We assume that the probability along the rapidity axis to 
find the hard gluon - fragmenting into a charmed quark 
pair - is defined by a triangle whose maximum is at ycm 
and which goes linearly to zero at y — asinh(xi-\/s/2TON) 
and y — ~asmh{x2\/s/2m^). The area of this triangle is 
set to unity so that it represents a proper probability. 

We have estimated the charmonium cross sections by 
the expression (taken from Ref. [23]) 



NN 



{s)^.Uail 



mi 



6{^s- 



'Soi 



(11) 



in which ^/s is the center-of-mass collision energy (per nu- 
cleon pair) and mi is the mass of the charmonium state i. 
The parameters a^lO, /3=0.775 and a are common for all 
charmonia and fitted to experimental data. The thresh- 
old factors read y^iol = rrii -\- 2m,N while the parame- 
ters fi were fitted separately in [23] for each of the states 
(/i=0.636, 0.581 and 0.21 for Xc J/V' and V'' j respec- 
tively). Above, the multiplicity label Xc denotes the sum 
of the three xoc Xic and X2c states. All these states decay 
into J/i/i and we have taken the sum of their branch- 
ing ratios (0.55) into i/ip as our branching ratio for 
the generic "Xc"- We have slightly re-adjusted the com- 
mon normalization factor a from 0.16 mb to 0.133 mb 
in order to reproduce exactly the total J/i/' production 
cross section a^Il — 3.0^6 measured by the PHENIX 
collaboration [46]. The rapidity distribution of i/%1) 's in 
proton-proton collisions at -^5=200 GeV, evaluated ac- 
cording to Eqs. (9), (10) and (11) is compared with the 
PHENIX measurement in Fig. 4. The agreement appears 
good enough so that we can estimate both the rapidity 
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density at mid-rapidity as well as the total cross section 
in our simple model. 



3.0x10" 



with the power-law function 



2.5x10" 



2.0x10"^ 
^ 1.5x10^ 



1.0x10" 



5.0x10" 




^ ^ (1 + [PT/Bfr, 



(12) 



-2-1 1 



see Figure 5. Here we have used the published data and 
fitted B=3.74 and n=5.11. 

According to our best knowledge, the transverse mo- 
mentum distribution of D mesons have not been measured 
in proton-proton collisions at RHIC beam energies. Pre- 
liminary data [48] in Cu + Cu collisions do exist as well 
as already published data [49,50] in d -|- Au collisions at 
■ysj^j^=200 GeV. The transverse momentum spectra of 
D mesons in these collision systems - divided by the num- 
ber of binary collisions - along the corresponding J/ip data 
inp+p collisions are shown in Fig. 5. The lines shown have 
the functional from of Eq. (12) and are fitted to the data. 
The D meson data do not allow to reliably fit both B and 
n (as well as the normalization) and so we have chosen 
to fix the parameter B = 3.74 as in the case of J/f/; and 
re- fitted n = 10.7 in order to describe the STAR data for 
d + Au collisions. 



Fig. 4. (Color on-line) Calculated and measured J/-!/) rapi- 
dity distribution in proton-proton collisions at ^/s—200 GeV. 
The data (full dots with error bars) are from the PHENIX 
collaboration [46]. The primary production of 3/ip, Xc and tp' 
are evaluated as described in the text. The "From decays" - 
line denotes the feeding from ip' and Xc states. The central 
rapidity region, in which the dilepton measurement is carried 
out, is indicated by the vertical lines. 



The reason we have taken the trouble to set up a model 
that can describe the J/"!/; rapidity distribution in p + p 
collisions is that we need this model to evaluate the ra- 
pidity distribution of D mesons. The open charm rapidity 
distribution is not yet measured at this beam energy and 
thus we need to calculate it. Since we have seen that our 
model can describe the J/ip data in this collision system, 
we can fairly safely assume that the same model will, at 
least approximately, describe the longitudinal part of the 
open charm momentum distribution as well. 

Let us turn our attention to the transverse directions 
now. As we discussed earlier, the transverse direction con- 
tributes most to the invariant mass of the dileptons from 
open charm decays at mid-rapidity. This is why we will 
rely on experimental data here. The transverse momen- 
tum distributions of D mesons are experimentally not well 
known, though. What is much better known is the rapid- 
ity and transverse mass distribution of J/ip 's in proton- 
proton collisions at RHIC. We deem that the momentum 
distributions of D mesons resemble the corresponding ones 
for J/-0 since the shape of the distribution - especially in 
the beam direction - should be primarily determined by 
the dynamics of the hadronizing charmed quarks. Thus, 
we assume that the D meson transverse momentum dis- 
tribution has a similar form as that for J/ip 's in the same 
collision system. The pr spectrum of J/ip [46,47] (mea- 
sured by the PHENIX collaboration) can be described well 
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Fig. 5. (Color on-line) D meson pT spectra per binary collision 
in d + Au [49] and Cu + Cu [48] collisions as well as the J/tp px 
spectrum [46] in p + p collisions. The lines are power-law fits 
to the data as described in the text. 



The dilepton continuum stemming from D meson de- 
cays attains an extra feature compared with the dilep- 
ton emission from other sources. Namely, since the elec- 
tron and positron are emitted by two different hadrons, 
the angular correlation as well as re-scattering effects can 
strongly alter the invariant mass spectrum of the final 
state dileptons. We assume here that the emitted leptons 
themselves always escape the collision zone unscathed and 
(re-)scattering effects can only take place on the hadronic 
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level. The angular correlations of the open charm hadrons 
are experimentally not well known and thus we will em- 
ploy theoretical estimates. 

We deem that the two charmed quarks are always emit- 
ted in 180 ° angle in their respective center-of-mass frame 
while this angle is typically much smaller in the labora- 
tory frame due to large Lorentz boosts in the longitudinal 
direction, especially at large forward and backward rapidi- 
ties. We assume here that the longitudinal direction is not 
different from the transverse directions in the CM frame 
of the fragmenting cc-pair, i.e. that Eq. (12) describes 
the joint distribution of any two momentum components 

"PT"^\/pi+Pl^ \/pI+pI=\JpI+pI in the CM frame 
of the charmed hadrons. The angular distribution among 
the produced D mesons is then taken exactly back-to-back 
correlated in their respective CM frame and we need to 
boost the momenta of the produced D mesons into the 
laboratory frame in order to evaluate the angular corre- 
lations in that particular frame. We have cross-checked 
our approach and verified that our angular correlations - 
evaluated as described above - agree well with correlations 
evaluated with the PYTHIA [51] event generator. Alterna- 
tively, we might also have adopted the angular correlations 
from the PYTHIA simulations. 

Wc have taken into account the 12 lightest D meson 
states {D^, D*^, D°, D^, D*^\ D*o, Df and D*^^) aU 
of whose mass is around 2 GeV. We have assumed that 
the relative primary multiplicity of these 12 states is de- 
termined purely by their spin -degeneracy while the total 
number of them is taken from the parametrization of Ref. 
[23]. On top of this, we have taken into account the em- 
pirical fact that in a jet fragmentation process, hadrons 
that include a strange or anti-strange quark suffer further 
suppression. We have used the canonical value 0.3 for such 
a strangeness suppression factor for the Dg and D* states 
in our analysis as in Ref. [23] which stems from PYTHIA 
calculations. One should not confuse this factor with the 
7s parameter included in the thermal model analysis. The 
factor 0.3 here concerns the hard scatterings only, while 
the 7s parameter takes into account also some of the soft 
physics on top of the suppression on the hard scattering 
level. 

The kinematics of the decays of the excited D mesons 
can be neglected because of the small mass difference be- 
tween the D* and D states. In each of the cases, the ex- 
cited state relaxes itself by emitting a very soft pion and 
thus the daughter D inherits practically the momentum 
of the decaying parent state. We deem that by properly 
considering all of the lightest D mesons (and their relative 
abundances) , we can then hope to extract the total charm 
production cross section in our analysis. Alternatively we 
could have just included the lowest lying D meson states, 
since these are the only ones decaying into leptons, but 
in this case the normalization factor, the number of final 
state D mesons would not have a clear physical interpre- 
tation as a cross section. Rather than taking the rapid- 
ity densities of different charmed hadrons at mid-rapidity 
as free parameters, we wish to estimate the total cross 
sections of the different states and suitably distribute the 



produced charmed hadrons at different rapidities. In order 
to estimate the rapidity density of D mesons around mid- 
rapidity, we take the total cross section for open charm 
production^ from the parametrization given in [23] and 
distribute these along rapidity as described earlier. 



5 Results for dileptons 

5.1 Proton-proton collisions 

Let us first look at the proton-proton collisions. Our cal- 
culated dielectron yields are compared with the experi- 
mental data in Fig. 6. The top panel of Figure 6 shows 
the low mass region where the emission is dominated by 
the decays of light mesons. One can see that the spectrum 
can be reproduced very well within the statistical had- 
ronization ansatz described in the previous sections. We 
mention that the p + p dilepton mass spectra are also well 
reproduced within the HSD transport approach [5] where 
the charm production and angular correlations have been 
evaluated within PYTHIA. 

The shapes of the two peaks arising from the direct 
decays of uj and (j) rnesons are essentially determined by 
the experimental mass resolution. The natural width of 
these peaks would be only about 8 and 4 MeV, respec- 
tively. The mass resolution of the experiments, however, 
exceeds the vacuum width of these hadrons and thus we 
have taken this into account by smearing the peaks accord- 
ing to Gaussian distributions with a width corresponding 
to an experimental mass resolution of 10 MeV for the light 
mesons except for the (f> -meson direct decay whose resolu- 
tion was taken to be 8.1 MeV [52]. The mass resolution of 
dielectrons from the charmed hadron decays was taken to 
be 2% of mg+e- . By doing so, especially the shape of the (j) 
meson peak is strongly modified and becomes compatible 
with the experimental data as seen in Fig. 6. We have not 
considered in-medium modifications of the spectral func- 
tions in this work and thus any apparent change in the 
vacuum spectral functions is solely due to experimental 
acceptance cuts. 

We find no discrepancy"^ between the data and our 
model in the low invariant mass region and can conclude 
that the data are very well described within the statisti- 
cal hadronization model assuming a thermalized fireball. 
The only non-equilibrium feature we have taken into ac- 
count here is the strangeness suppression factor 75. By 
this factor we introduce a new - partly free but corre- 
lated - normalization for hadrons consisting of one or more 



^ acc—o-j;,+ +a]jo+ o-D^+cr£,* + +crj5.o+CTDj=485 fib 
'^ The invariant mass region between the co and 4> meson 
peaks is somewhat under-estimated in our approach. We deem 
that this is due to our approximation that the angular corre- 
lations among the D mesons are back-to-back in the center-of- 
mass frame. In reality, a backward peaked but narrow angular 
distribution is expected, which softens the spectrum somewhat 
and populates the low invariant mass region. This effect, how- 
ever, would be only visible in the region between the cu and (jf> 
mesons. 
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Fig. 6. (Color on-line) Invariant mass spectrum of pairs of 
electrons and positrons in proton-proton collisions at -^s = 
200 GeV. The data are from the PHENIX collaboration [1] 
while the contribution from different dilepton emitting sources 
are calculated as explained in the text. The full thick black line 
denotes the sum from all relevant sources. The LMR is shown 
in the upper panel while the whole invariant mass range is 
shown in the bottom panel. 



over-estimated. The 7s parameter affects also the yields 
of T] and 77' mesons, whose primary thermal production 
rates are multiplied by 7^ in order to take into account 
the fact that these mesons are considered to carry "hid- 
den strangeness" . Thus the value of 7s is mostly (but not 
solely) determined by the </> meson. We remind the reader 
here that, besides the overall normalization, we have not 
fine tuned any of the thermal parameters in this work, in- 
stead, we have used the values fitted in Ref. [27] to STAR 
data. 

The whole invariant mass spectrum is shown in the 
lower panel of Figure 6. Let us look at the IMR between 
the (f> and J/ip peaks now. The solid and dashed black 
lines are evaluated as explained in the last Section im- 
plementing the transverse momentum profile fitted to the 
D meson transverse momentum spectrum in d -I- Au colli- 
sions. The only difTerence between these two lines is that 
the upper solid line is evaluated assuming back-to-back 
angular correlations in the center-of-mass frame between 
the two fragmenting D mesons while the lower dashed line 
is the same with random correlations. From the figure it 
is clear that the model with random correlations does not 
describe the IMR spectrum properly for p + p reactions. 

From this we can conclude, in accordance with the 
original PHENIX publication [1] , that a model with strong 
correlations (no final state interactions) among the pro- 
duced D mesons seems to be favored by the data over 
the random correlation in case of proton-proton collisions. 
We mention that in the extreme case of exact back-to-back 
correlations in the laboratory frame among the D mesons, 
the slope of the dilepton continuum in the IMR is repro- 
duced but in this case too few dileptons are emitted in 
the LMR. We do not explicitly show these results but 
note that the continuum between the w and would be 
under-estimated by about an order of magnitude in this 
over-simplified scenario. 

We may conclude that once all relevant kinematical 
as well as acceptance effects are taken into account, the 
electron-positron invariant mass spectrum can be under- 
stood very well in both the LMR and IMR in proton- 
proton collisions within our simple model and parametriza- 
tion. This will serve as a baseline for our comparative anal- 
ysis in heavy-ion collisions for different centralities. We 
mention in passing that the dilepton mass spectra from 
our simplified model agree with those from the HSD ap- 
proach [5] on the 20% level when using vacuum spectral 
functions for the hadrons in the transport approach. 



5.2 Heavy-ion collisions 



strange quarks. In practice, the action of the 75 parameter 
is best visible in the cj) meson peak. Since there are no reso- 
nances decaying into the (j) meson the sole rapidity density 
of the meson is calculated according to Eq. (1). With 
our choice for the value of 7s the meson rapidity den- 
sity is multiplied by a factor 7|.=0.36. Obviously, without 
this extra strangeness suppression the production of dilep- 
tons from the decays of (jj mesons would be dramatically 



Let us turn to heavy-ion collisions now. We have seen in 
the previous Section that the dilepton invariant mass spec- 
trum can be well understood in proton-proton collisions 
and we will now try to extrapolate our results from p + p 
collisions to heavy-ion collisions to identify the magnitude 
of possible additional production channels from partonic 
sources as suggested e.g. in Ref. [53] for SPS energies. 

The non-charmed hadron yields are expected to scale 
with the participant number (A^part) when comparing p-l-p 
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and heavy-ion collisions while - assuming that the charmed 
quark pairs are created solely in the initial hard scatterings 
- the number of binary collisions (A^bin) should be the 
correct scaling variable for the charm sector. Both TVpart 
and iVbin can be estimated within the Glauber model and 
are conventionally used also in the experimental analyses. 
We have used the same values for A'part and iVbin [31,54] 
that are employed by the experiments at RHIC in order 
to be consistent with similar previous analyses. 



5.2.1 Low invariant mass region 

The Au + Au and Cu + Cu collisions are treated different 
in this paper. For the Au-\- Au collisions we have fitted the 
width of the clusters' transverse momentum distribution 
and the system volume to the PHENIX data, see Fig- 
ures 2 and 3, and evaluated the dilepton invariant mass 
spectrum in Au + Au collisions implementing the 'opti- 
mized' parameters. Detailed measurements of identified 
hadron transverse momentum spectra are not yet avail- 
able in the Cu + Cu collisions and thus we have evaluated 
the yields of dielectrons in the LMR in Cu -\- Cu collisions 
such that we use the transverse momentum profile fitted 
to the Au -\- Au collisions also in Cu + Cu collisions at the 
same iVpart- The volume in Cu + Cu collisions at different 
centralities is evaluated by scaling the system volume in 
p + p collisions with the number of participants. As one 
can see from the Figure 3, the fitted volumes scale with 
A'part in the Au + Au collisions in the A'part region relevant 
for Cu + Cu collisions and thus we deem that the iVpart 
scaling should be a good approximation in the Cu + Cu 
case, at least in the non-central collisions. 

We have evaluated the dilepton yields in Au + Au and 
Cu + Cu collisions as described above, see Figs. 7 and 8. 
In the lower panel of Fig. 7 we show our results together 
with the experimental data [55] in minimum bias Au + 
Au collisions in the whole invariant mass range while in 
the top panel a zoom to the LMR is presented. In Fig. 
8 the experimental results from Refs. [2,56] in different 
centrality bins both in Au + Au (top panel) and Cu + Cu 
(bottom panel) collisions are shown. The Cu + Cu data in 
different centrality bins are still preliminary and the error 
bars are not yet available. 

In the following we will concentrate on discussing the 
Au + Au collisions. Essentially the same conclusions will, 
however, hold also for the Cu + Cu collisions. Let us take a 
closer look at the LMR first. From the top panel of Fig. 8 
one can see that the data can be described well in the most 
peripheral centrality bin in Au + Au collisions in the LMR. 
This centrality class is special among the centrality classes 
in Au + Au collisions because in all other centrality classes 
the relative strangeness production is found to be nearly in 
chemical equilibrium with 75 w 1 [19]. We have taken the 
increase in relative strangeness production into account 
in our analysis by setting the 7s parameter to unity in 
every other centrality bins except in the most peripheral 
collisions, i.e. in the centrality bin labeled with "5" in 
Fig. 8a and "4" in Fig. 8b. The resuhs with both 75 = 
0.6 (solid) as well as 75 = 1 (dashed) are shown for the 
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Fig. 7. (Color on-line) Invariant mass spectrum of pairs of 
electrons and positrons in minimum bias Au + Au collisions at 
v'iNN=200 GeV [55] compared with model calculations. The 
low invariant mass region is shown in the top panel while the 
results in the full invariant mass region are shown in the bot- 
tom panel. The solid line indicates the model results scaled 
from p + p collisions such that the charmed hadron yields are 
scaled with the number of binary collisions A''bin while the non- 
charmed hadron yields are scaled with 110 1/pp (the volume in 
p + p collisions). The dashed curve (bottom panel) indicates 
the results scaled from p + p collisions assuming a random cor- 
relation among the D mesons. The calculation with the 75 
parameter set to unity is indicated by the dashed line (top 
panel) slightly above the solid one in the LMR. The dotted 
curve (bottom panel) indicates the results evaluated with an- 
gular correlations kept the same as in p+p collisions. The solid 
and dotted error bars stand for the statistical and statistical 
-I- systematic errors added in quadrature, respectively. 
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Fig. 8. (Color on-line) Invariant mass spectrum of pairs of 
electrons and positrons in Au + Au (top) [2] and Cu-\-Cu (bot- 
tom) [56] collisions at y^f^j^=200 GeV in different centrality 
classes compared with the model calculations. The centrality 
bins are labeled from central to peripheral as 1 (0-10%), 2 (10- 
20%), 3 (20-40%). Centrality bin 4 consist on (40-94%) and 
(40-60%) most central events in Cu + Cu and Au + Au colli- 
sions, respectively, while the centrality bin 5 includes (60-92%) 
most central collisions. Both the data and model calculational 
results are scaled with factors of 10'' (0-10%), 10^ (10-20%), 
1 (20-40%), 10"^ (40-60%) as well as 10"^ (60-92%) and (40- 
94%) for clarity. The double dotted lines indicate the model re- 
sults scaled from p + p collisions such that the charmed hadron 
yields are scaled with A''bin. The non-charmed hadron yields 
are scaled with A^part in the case of Cu + Cu collisions while 
in Au + Au collisions we have scaled the yields with the fit- 
ted volumes. The dashed lines indicate the results scaled from 
p+p collisions with random correlations among the open charm 
hadrons. The starting point of the downward pointing arrows 
denote the data points which are defined as upper limits only. 



niinimum bias Au+Au collisions in the top panel of Figure 
7. One can see that the increase in relative strangeness 
production as a function of centrality is not strong enough 
to explain the excess in the LMR in minimum bias Au+Au 
collisions. 

From Figure 8, one can see that the LMR data in 
Au + Au collisions in peripheral and in semi-central bins 
can be reasonably well described within the statistical 
hadronization model. On the other hand, in the two most 
central bins as well as in the minimum bias collisions the 
increase in strangeness production can not explain the ex- 
cess of dileptons in the low invariant mass region from 0.15 
to 0.6 GeV. In Cu + Cu collisions, it seems that there is 
significant excess in the LMR over the hadronic cocktail 
only in the most central collision bin while the LMR is 
fairly well described for the other centralities. 

We have, furthermore, studied the effect of the trans- 
verse flow on the dielectron invariant mass spectrum in 
Au+Au collisions by comparing the results evaluated with 
the maximum and minimum width for the clusters' trans- 
verse momentum distribution. The largest (see Fig. 2) 
width, crpy = 13.5 GeV, was fitted in (15-20%) most central 
collisions while for the p + p and (60-92%) most central 
collisions we use crpj,=10.0 GeV. We have evaluated the 
dielectron invariant mass spectrum with these two widths 
keeping all their parameters fixed (T = 170 MeV, V = Vpp 
and 75 = 1) and calculated the ratio of the resulting in- 
variant mass spectra of dielectrons in the PHENIX accep- 
tance, see Figure 9. 
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Fig. 9. (Color on-line) Ratio of invariant mass spectrum of 
electron-positron pairs within the PHENIX acceptance evalu- 
ated with two different transverse flow profiles keeping other 
parameters fixed. 



As expected, the increase in the transverse flow en- 
hances the dielectron yields in the PHENIX acceptance 
especially at very low invariant masses while above the w 
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meson peak equal amounts of dileptons hit the PHENIX 
acceptance with both transverse flow profiles. In general 
the effect is moderate and the increase in flow can increase 
the dielectron yields up to 30% in the 7r° region, much less 
than the observed excess by factors up to 4-5. Thus, one 
can conclude that the broadening of transverse momen- 
tum distributions as a function of centrality plays only a 
minor role in the dielectron radiation from the decays of 
light hadrons. 



5.2.2 Further correlated sources in the LMR 

Besides the decays of light mesons, there are other corre- 
lated sources for dielectron production in high-energy nu- 
clear collisions. In this Section, we will address some but 
not all of such processes. The contribution from most of 
these channels is very small and thus we concentrate here 
on the dominant channels, only. These additional chan- 
nels, similarly to the the D meson case, arise in corre- 
lated decays of light hadrons that do not directly decay 
into electrons but produce single electrons via intermedi- 
ate hadrons. An example of such a process is the (/)-meson 
decay — > K^K~ . The (j) meson has a long lifetime and 
a small hadronic cross section and thus many of the (j) 
meson decays take place outside the fireball. Now it can 
happen that the kaons from the (j) meson undergo semilep- 
tonic decays like K^ — >■ n'^e^Ve and charge conjugate for 
the K~ . If the (j) meson decays outside the fireball, the 
correlations are preserved and the experiment measures 
essentially additional correlated dielectron radiation from 
the (j) meson which survives the experimental subtraction 
procedure from uncorrelated e^e~ decays. 

Unlike the (p meson, most of the short living resonances 
are expected to decay inside the fireball thus destroying 
the correlated signal. However, some of the interactions 
do take place in the dilute corona of the fireball in which 
case the correlated signal can be preserved even in central 
heavy-ion collisions. We do not perform a precise calcula- 
tion for these correlations and make an estimate for the 
relative magnitude of core and corona emission, instead. 
Our results thus have to be considered as an upper limit 
estimate for a correlated background emission. 

We have considered this type of correlated dielectron 
radiation from a selection of fairly light strange and neu- 
tral hadrons whose branching fractions into KK are size- 
able and reasonably well known. The additional channels 
we have studied in this work are listed in Table 2. We con- 
sider the kaonic channels only in this work, even though 
most of the processes could proceed via the tttt channel 
as well. The pionic channels are found to be sub-leading 
compared to the kaonic channels and thus we have omitted 
pionic channels in this work. 

We have considered all K'^'s as fixed 50% - 50% mix- 
tures of K^ and K^ states and ignored the time dependent 
neutral kaon oscillations. Making a distinction between 
the long and short living states is important since the 
K'j^ has roughly a factor of 100 larger probability to de- 
cay semileptonically than the K'g. Electrons and positrons 
stemming from X^ decays are taken into account only if 



Hadron 








/o(980) 
/i(1285) 
/2(1270) 


K+K- 
K+K- 




KK-K 


/o(1350) 
/i(1420) 
/2(1525) 


K+K- 

K*+K- + c.c 

K+K- 


K*°K'' + c.c 


- 


/o(1500) 
/i(1510) 
/2(1430) 


K+K- 

K*+K- + c.c 
K+K- 


K'-'RO 
K'^K'' + c.c 


- 


a!j(980) 
K{m2f 


K+R- 
K+K- 


K'-'RO 


- 



Table 2. List of additional decay channels relevant for dielec- 
tron production in high-energy collision systems. 



the K^ has decayed before the first detector (2 meters 
from the primary vertex) which reduces significantly the 
di-clcctron yields from the K'^ decays. 

The contributions from the additional correlated chan- 
nels in (0-10%) most central Au + Au collisions are shown 
in top panel of Figure 10. We have added some of the 
channels in Figure 10 for clarity. The K* denotes the sum 
of K*^, K*^ and K*^ while /012 is the sum of /o, /i and 
/2 and similarly /g^2 denotes the sum of the /o(1370), 
/i(1420) and /2(1525). The contribution from the /o(1500) 
is small and is not shown. For comparison, also the hadronic 
cocktail contribution from 77 and p'^ mesons are shown by 
the thick solid lines. One can see that the correlated chan- 
nels might indeed give a sizeable contribution to the di- 
electron invariant mass spectrum in the low invariant mass 
region and might even over-shine the standard hadronic 
cocktail emission precisely in the invariant mass region 
where the large excess was measured by PHENIX in cen- 
tral nucleus-nucleus collisions. 

In the bottom panel of Figure 10 the standard hadronic 
cocktail result (solid) and the cocktail + additional corre- 
lated emission (dashed) are compared with the PHENIX 
data in the most central bin. Our results show that the 
correlated background from the exotic mesonic states may 
result in a large enhancement of the dielectron yields in 
the LMR, although re-scattering effects could significantly 
alter the results at least for some of the correlated chan- 
nels. The sizeable contribution - relative to proton-proton 
reactions - is due to a factor 7I which increases roughly 
by a factor of three from peripheral to central nucleus- 
nucleus collisions. Note, however, that we have addressed 
upper limits, only, and that a decorrelation by interactions 
in the hot medium is expected. Nevertheless, even our up- 
per limit is clearly below the PHENIX signal for central 
nucleus-nucleus collisions and we may conclude that the 
additional channels considered here should not be respon- 
sible for the dilepton excess seen experimentally. 
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Fig. 10. (Color on-line) The upper limits for the invariant 
mass spectrum of pairs of electrons and positrons in (0-10%) 
most central Au -\- Au collisions arising from various correlated 
sources (top panel) in comparison with the PHENIX data and 
the standard cocktail calculation (bottom panel). 



5.2.3 Intermediate mass region 

We have estimated the charmed hadron yields in heavy- 
ion coUisions by multiplying the corresponding yields in 
proton-proton collisions with the number of binary colli- 
sions for the contributions from charm mesons. This pro- 
cedure should hold for open charm hadrons but actually 
we know that the production of charmonium states suffer 
suppression in the hot and dense environment of partonic 
and hadronic nature. This suppression is usually expressed 
by the ratio Raa = jv dk^^-^ which parametrizes the de- 
viation from the simple A'bin scaling. PHENIX has mea- 



sured this quantity for J/-0 's both in Au + Au [57] as well 
as in Cu + Cu [58] collisions. So far there is no fully con- 
vincing model calculation that explains the observed Raa 
on a satisfactory level and thus we have taken the simplest 
approach and assumed that Raa is a linear function^ of 
centrality. This approximation seems to hold sufficiently 
far away from the ends of the whole centrality range. For 
alternative curves for the Raa with centrality we refer the 
reader to the review [23]. 

The different model results in Figs. 7 and 8 take into 
account the charmonium suppression. It is clear from Fig. 
8 that it is necessary and sufficient to include the charmo- 
nium suppression effects in order to describe the i /ip peak 
correctly at all centralities. Since the Raa ratios of the 
excited charmonium states are not yet measured, we have 
not implemented any correction for the ip' . 

Let us now turn our attention to the slope between 
the (j) and i /ip peaks. From the discussion above it is clear 
that both the very low invariant mass (tt*^ Dalitz decay) 
as well as high invariant mass {J /^) regions are fairly well 
described by the model. This is a good starting point to 
address the physics of the IMR in between. 

In p + p collisions there is no medium that would 
distort the correlation among the emitted open charm 
mesons and indeed, we have seen that the p + p data 
can be best described by assuming strong correlations 
among the D mesons. In central heavy-ion collisions on 
the other hand, there are several hundreds of hadrons 
emitted in each event and one would expect that the 
produced charmed hadrons interact with the surrounding 
medium thus destroying the initial correlations. We have 
studied the reinteractions of the charm mesons quantita- 
tively within the HSD transport approach (cf. Ref. [23]) 
and have calculated the probability that neither of the two 
D mesons interacts with the surrounding medium. In this 
and only in this case the angular correlations among the 
emitted electrons and positrons would be preserved and 
remain similar to the p + p case. 

Our results from the HSD calculations are shown in 
Fig. 11. The probability that the angular correlations re- 
main the same a.s in p+p collisions is calculated as a func- 
tion of collision centrality in Cu -\- Cu (open circles) and 
in Au + Au (filled circles) collisions at -ysj^j^=200 GeV. 
For practical purposes, we have parametrized these prob- 
abilities by the functions shown in the Figure. In the case 
of minimum bias Au + Au collisions, a proper weighted 
average over the whole centrality range (denoted by the 
filled square in Fig. 11) is used in our calculation instead 
of the explicit parametrization. 

We can now study the invariant mass spectrum of the 
dilepton continuum in heavy-ion collisions in a more re- 
alistic scenario in which the slope of the IMR arises as a 
superposition of correlated and un-correlated open charm 
decays. All the three cases are shown for the minimum bias 
Au + Au collisions in the lower panel of Fig. 7. The result 
retaining the correlations like in p + p collisions is shown 
by the double dotted line in Fig. 7b, while the random 



* i?4 4-*" = 0.75 - 0.00137iVp 
Raa 






= 1.0-0.00516Afpart 
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Fig. 11. (Color on-line) Probability that neither of the 
D and D mesons collide with the surrounding medium in 
Cu -I- Cu (open circles) and Au -\- Au (closed circles) collisions 
at v^NN~200 GeV as evaluated with the HSD transport ap- 
proach. The filled square denotes the value for minimum bias 
Au-\- Au collisions. The lines are parametrization that we have 
implemented in our calculations. 



In addition to the change in the angular correlations 
among the open charm mesons, also the magnitude of the 
(transverse) momentum of the D mesons can change due 
to interactions in the fireball. The PHENIX collaboration 
has measured the Raa of single electrons coming from 
heavy quark decays [64] in Au -I- Au collisions and found 
that this Raa is compatible with unity up to pt of 2 
GeV at all centralities. Thus, heavy-quark energy loss can 
(significantly) modify the di-clcctron yields in the IMR 
coming from D meson decays only at transverse momenta 
larger than 2 GeV. 

In order to estimate how this affects our results, we 
have calculated the single electron -f positron transverse 
momentum spectrum within the PHENIX acceptance and 
divided the results in four classes according to the corre- 
sponding invariant mass of the di-electron pair Mg+g-. 
The pt spectra are shown in mass windows of Mg-i-g- € 
[0,1] ; [1,2] ; [2,3] and [3,4] GeV in Figure 12. The solid 
lines indicate the results in the fully correlated case while 
the dashed lines are evaluated in the random correlation 
picture. All curves are normalized to unity. One can see 
that a significant fraction of electrons and positrons with 
Pt >2 GeV contributes to the invariant mass spectrum at 
invariant masses larger than 3 GeV, only. Our calculations 
concentrate on the region Mg+g- G [0,Mj/^] and thus we 
can conclude that our results are only little affected by the 
heavy-quark energy-loss effects in the medium. In partic- 
ular, the invariant mass region Mg+g- S [1,2] GeV - in 
which our calculations under-estimate the measurements 
- appear little affected by the heavy-quark energy loss. 



correlation case is represented by a dashed line. The solid 
line in between shows the invariant mass spectrum assum- 
ing that in 33% of the cases the dilepton pair stemming 
from the open charm decays retains the initial correlations 
while in 67% of the cases at least one of the D mesons has 
scattered and thus the correlation is destroyed. One can 
see that the most realistic case naturally interpolates be- 
tween the two extreme cases and seems to describe the 
experimental data best (within the error bars). 

Some further information on the IMR is gained by hav- 
ing an explicit look at the centrality dependence of the 
dilepton yield for Au + Au and Cu + Cu collisions at the 
top RHIC energy. The solid lines in Figure 13 are eval- 
uated in this mixed scenario - based on the HSD rescat- 
tering calculations - and one can see that this approach 
practically underestimates the experimental spectra in the 
IMR for both systems. 

Presently, we may only speculate that there seem to be 
further channels of possibly partonic nature in the IMR as 
suggested by several groups independently. Within ther- 
mal models this excess might be addressed as thermal 
dilepton radiation from the QGP [59] while also hadronic 
p+p or 7r-|-ai scattering might contribute as suggested by 
van Hees and Rapp [7,60]. The studies by one of the au- 
thors on explicit partonic reaction channels in Refs. [61, 
62] allow for an implementation in the PHSD transport 
approach [63] which hopefully might clarify this issue in 
the near future. 



6 Discussion 

As we have seen in the previous Section a scaling of the 
dielectron yields from p + p to heavy-ion collisions can 
lead to a surprisingly good description of the data in the 
peripheral collision systems but fails for the more central 
collisions. The question we wish to address here is that 
if we can, nevertheless, understand the observed excess 
in terms of hadronic degrees of freedom or if additional 
partonic productions channels have to be incorporated. 

According to the statistical hadronization model fits 
to p + p and Au + Au collisions at RHIC, the intensive 
thermal characteristics of these systems seem very simi- 
lar at mid-rapidity. All of the light mesons decaying into 
dileptons are completely neutral and so the dilepton pro- 
duction rate does not actually depend on the chemical 
potentials at all and thus, besides the 75 parameter dis- 
cussed before, the temperature is the only intensive pa- 
rameter left in the model that could lead to a non-trivial 
scaling behavior seen in the data. However, one of the 
lessons we have learned from the SHM fits to RHIC data 
is that the temperature is the same in heavy-ion collisions 
at all centralities and this temperature coincides with the 
one extracted from the p+p collision data. We have, never- 
theless, checked that a moderate increase (T = 170 — )■ 190 
MeV) in temperature can not explain the observed excess 
in the LMR. One should notice that due to the momen- 
tum cut pt > 0.2 GeV the increase in temperature af- 
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Fig. 12. (Color on-line) Transverse momentum distribution 
of electrons and positrons (within the PHENIX acceptance, 
normalized to unity) coming from semi-leptonic decays of 
D mesons in four different invariant mass windows M^+^-. 
The solid lines indicate the spectra in the fully correlated case 
while the dashed lines arise in the random correlation picture. 
Electrons and positrons with pr larger than 2 GeV (right from 
the solid vertical line) exhibit energy loss while propagating 
through the hot and dense medium [64]. 



fects more prominently the dilepton yields from the vector 
mesons than the emission from tt" and thus a change in 
temperature affects different regions in the invariant mass 
spectrum with different strength. One would nevertheless 
need unreaUstically large temperatures of T > 200 MeV 
if one attempts to assign the dilepton excess (seen by the 
PHENIX Collaboration) to an increase of the fireball tem- 
perature. We rule out such a possibility. 

Independent previous model calculations [5,?, 7] have 
been compared with the PHENIX data in the original pub- 
lication [2] and we refer the reader to Figs. 41 and 42 of 
Ref. [2] for details. Our results agree qualitatively and 
also quantitatively with the previous model calculations 
in that the proton-proton collisions are well described in 
the whole invariant mass range while none of the analyses 
can explain the excess in the central heavy-ion collisions in 
the low invariant mass region. In extension of the previous 
studies we have, furthermore, investigated the possibility 
that the observed excess might stem from further semilep- 
tonic correlated kaon decays which are enhanced in central 
nucleus-nucleus collisions relative to proton-proton reac- 
tions by roughly a factor of three. Our upper limits for the 
dominant channels considered here clearly show that also 
these additional 'background sources' are not responsible 
for the large excess seen by the PHENIX Collaboration in 
central heavy-ion reactions. 

The charmed sector or the IMR (if considered) has 
been treated essentially in a similar fashion by the PHENIX 
Collaboration as in this work and no solid conclusions have 
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Fig. 13. (Color on-line) Invariant mass spectrum of pairs of 
electrons and positrons in Au-\-Au (top) [2] and Cu + Cu (bot- 
tom) [56] collisions at y^j^j,j=200 GeV in different centrality 
classes compared with model calculations. The centrality bins 
are labeled from central to peripheral as 1 (0-10%), 2 (10-20%), 
3 (20-40%). Centrality bin 4 consist of (40-94%) and (40-60%) 
most central events in Cu + Cu and Au + Au collisions, respec- 
tively, while the centrality bin 5 includes (60-92%) most cen- 
tral collisions. Both the data and model calculational results 
are scaled with factors of 10* (0-10%), 10^ (10-20%), 1 (20- 
40%), 10"^ (40-60%) as weU as 10"^ (60-92%) and (40-94%) 
for clarity. The solid lines indicate the model results scaled from 
p + p collisions such that the charmed hadron yields are scaled 
with A^bin while the non-charmed hadron yields are scaled with 
A'^Part (Cu + Cu) and with the fitted volume [Au + Au) while 
correlations among the open charm hadrons are evaluated in 
the mixed procedure as described in the text. The starting 
point of the downward pointing arrows denote the data points 
which are defined as upper limits only. 



16 



Manninen et al.: Dilopton production in p + p, Cu + Cu and Au + Au collisions at ^/s = 200 GeV 



been possible, so far. In order to go beyond the previous 
attempts we have calculated the rescattering probabili- 
ties of charm mesons dynamically (within HSD) which 
allows to estimate the amount of uncorrelated electron + 
positron pairs from D meson decays as a function of the 
centrality of the reaction. Our final results for Au + Au 
and Cu -f Cu suggest that we clearly underestimate the 
preliminary yield from PHENIX which might point to- 
wards partonic sources - as suggested in Ref. [62] - in the 
intermediate mass regime. 

So far we have considered only events with exactly one 
charmed quark pair. Processes leading to un-even amounts 
of charmed quarks are possible but they are more rare 
than the case we have studied and the corrections are 
probably not very large. Events with, e.g. 3 D mesons 
tend to populate the low invariant mass region, in which 
the open charm contribution is insignificant, because in 
that case only one of the possible two "dilepton" -pairs is 
(strongly) correlated. To finally clear up the situation we 
are going to carry out non-perturbative calculations on 
correlated charm dynamics within the PHSD transport 
approach [63] that also includes the dynamics of charm 
quarks in the partonic phase. 
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